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Abstract 



The FORG opeeration principle uses a recirculating ring resonant cavity to get a rotation- 
induced Sagnac effect enhancement [1]. It grants to a FORG a comparable sensitivity in relation 
an I-FOG 2 that has the fiber lenght 9/2 times longer, where 9 is equal its finesse factor. Other 
advantages is despite of thermal drift because the FORG uses less quantity of fibre than the I-FOG, 
giving to the first less thermal drift sensitivity than the last. But, due to the Kerr effect and couplers 
loss, that causes parasitics effects like dissimetries on signal response [2] and cross polarization, due 
to these facts, the FORG complexity is increased a lot. The signal response dissimetry caused by 
kerr effect can be corrected by phase nulling method, was proposed by Hotate [3]. The proposal of 
this work is to show a applied to a FORG technique that simplifies the signal treatment, employing 
all digital setup, like as filter banks and wavelets methods, resulting in a maximally flat scale factor. 
In this investigation are presented over the simulations results, employing the modified digital FM 
spectrometry techniques by decimation and interpolation techniques over a ring resonator that pursuit 
a 10 meters SM-PM lenght fiber coil and 10 centimeters of diameter, with a 1.55/um laser source. The 
advantages of these techniques are to simplify the electronic circuitry, offering an upgrade facility, 
using only one DSP (Digital Signal Processor), realizing all needed functions. The investigation of 
this method is based in a optical field switching scheme and digital frequency domain spectrometry. 
The purpose of this work is to describe this digital technique, well as the simulation results, discussing 
about this technique use and its limitations. 



1. Introduction 



The FORG operation principle is based on a recirculating light into the resonant cavity or a fibre 
ring resonator. It increases its sensitivity, enhancing the rotation induced Sagnac effect [1] by a 9/2, 
where 9 is equal its finesse factor. In this work we intend to show one easy solution to signal process 
employing less hardware. This technique employs light switching and digital frequency decimation 
methods to get the rate measurement signal. 



2. Theoretical Background 



The proposed setup is depicted in figure 1, that uses two electro-optic modulators, one in each 
Mach-Zehnder interferometer branch, each of wich working as an amplitude modulator (see ref.[9] 
for more details). Jointly of these modulators, have two electronic switches working syncrinizatelly. 
The polarization controllers and optical isolators had been suprcssed of fig.l for simplicity. The set 
of equations that describe the FORG theory will be given in the following pharagraphs. The FORG 
response equations can be found in reference [4], which will be given below. 
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Where, in eq.2, AL' R (Q) = L ^ n , and oq = j^. The values 70, are, respectively, the factional 
coupler intensity loss and resonant coupling coefficient (re r = (1 — 7o)~ 2a ° L , where «o is the exponential 
attenuation per lenght unit). The index % in Ri((3) denotes the i th port (i = 1,2) output. These 
equations are derived from coupling matrix, depicted in equation 5: 



1 Fibre Optic Resonant Gyroscope 

2 Tm.tp.rfp.mm.pt.rir. Fihrp Onti.r dumsmnp 



CXJ 

A{C + AB 2 a [aAC] m " 1 e- J '(- woroT ^r=o*"(*- fcT ))} 



m=l 

oo 



En 



A{C + AB 2 a MC] m_1 e- J ' ( - WomT ^r=o 



0+(t-fcr)) 



} 



(3) 



(4) 



m=l 



Where the A, B, C and D constants are the complex elements of the unitary coupling matrix, and 
a is equal the attenuation due to the fiber lenght (a = exp[— q.qL\). The coupling matrix is shown 
below 



C = (1 - 7o) : 
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C D 
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The A and B coefficients are equal to (1 — /e)5, and C and D coefficients are equal to j/ca. The 
kappa constant is the intensity coupling coefficient. Taking account the phase modulator character- 
istics (e.g.: an LiNbOs electro-optic modulator or a fiber wrapped and bounded in a piezoelectric 
cylinder), its response has the following equation: 
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Changing the term in summation into the exponential in eq. 3 and 4, by eq. 6, after some algebraic 
manipulations, these expressions can be given by (see reference [7] for more details): 
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To get a maximum sensitivity rotation response(and simplifying the computational cost), the 
equations 7 and 8 must obeys the following condition: 
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The resulting frequency in eq. 9 is the necessary frequency to reach the maximum phase variation. 
Expanding the expression 2 in Taylor's series, into the discrete-time domain (t — > nT s ), changing the 
values into the sine, already depicted in eq. 1, taking account the frequency value in eq. 9, we get 
the following expressions: 
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Where the Tj,, V 2 , and T\ coefficients are, respectivelly: 
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The 4> function is equal |8 m sm[^], where (3 m is the modulation index, and is the discrete 
modulation frequency. The A', B', C and E' Coefficients of 12 to 15 are described as 



A' 
B' 

c 

E' 



l-7o 
(1 

(1 + K0 

4ft r 



(16) 
(17) 
(18) 
(19) 



The sine and cosine terms in eqs. 12, 13, 14 and 15 are expandable in I s * kind of Bessel's series. 
These series, in sine and cosine expansion, are given by: 
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The arguments into the summation symbols in equations 21, 23 and 25 and 26 are digital 
modulation frequency W m multiple integer, where W m = where N = j 2 -. These arguments are 
made equal a multiple integer of ^ by decimation [8] to reduce the output ripple factor signals due 
to several generated harmonics by nonlinearity signal response during the phase modulation. The 
decination factor best choice is made equal to ^. Then, the 21, 23 and 24 become: 
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The equations 27, 28, 29 and 30, after some algebraic manipulations, can be simplified to the 
following equations (see reference [7]): 
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The sin(m/3 m ) and cos{m(3 m ) values (m integer) are ever constants. Making AP(f2) equal 
P CCTO (J7) — P CW (Q), and replacing the terms of <j>(k) sines and cosines of equations 12 to 15 by equa- 
tions 31 to 34, the even terms are eliminated, giving a signal nonreciprocity, i.e., the rotation sense 
information. The resulting equation is given by: 

AP(Q) = 2[Tl(7T)AL' R (n)a ) + r|(7r)AL' R (^ ) 3 + 0(5)] (35) 

The measured value of Q is recursively computed using the priori values (at the k — 1 instant). 
Then, we have the priori and posteriori expressions, resulting the following expressions: 
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Changing to the matricial form and neglecting the higher order and considering that the data 
aquisition time lack is very short to consider the among of f2 variation (the processing time is shorter 
than the variation), the equations 36 and 37 can be put in the matricial form 
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To computer the measured value of Q, after some matrix manipulations, the equation 38 can 
be written in the following equation 
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The modulator factor choice is very important because the unsuitable (3 m value can cause a great 
overshoot, turning the system critically damped. In this work, the (3 m is determinated by simulation. 
To correct the modulation factor, is needed to adjust the phase modulation driving amplifier gain. 
Note that the output matrix provides proportional to Q k and Qf signals, where the Q k quantity is 



3. The Setup and Operation Principle 



The operation principle is based on syncronization field intensity combination, where the coun- 
terpropagating fields intensity are digitally processed. These signals are sent to a processing block 
separatelly into the time domain by the optical switch and electronic switch. Then, these signals are 
combined by difference between them, cancelling the even powers of Qr. The frequency decimation of 
these signal have the finallity of eliminate the modulation frequency harmonics that simplifies the pro- 
cessing alghoritm. The system setup is depicted in figure 2, where : PD=photodiode, LD=laser diode, 
PM=phase modulator, OS=optical switch, OC=optical coupler, BS=beam splitter, FOC=fibre-optic 
coil, PC=polarisation controller, OI=Optical isolator. 

The correlation effect and the offset can be corrected by an additional photodiode elimination 
and the employing of an optical switch pair (shown in figure 3), working syncronized with two CMOS 
switches together by a clock pulse train (the optical switch operation principle is described in reference 
[9]). These optical switches can be integrated in a LiNbO% substrate. Each optical switch allows each 
field pass in a time interval equal of -y to the photodiode. The electrical switch scheme and electrical 
signal conditioning setup is proposed in figure 4 and 5, respectively, where Ai^ is the generated signal 
difference by counterpropagating fields. 

These optical switches provides a syncronizated signal difference with two CMOS switches, that 
separate the I c w and I c cw from photodiode at several instants of . Note that the I c cw is delayed by 
a time lack of ^f. It is necessary that the I c w and I c cw signals arrive at same instant to the subtractor 
to avoid the signal nulling in equation 35. These signals is subtacted one of each and decimated in 
frequency, explaned already in equations 27 to 34. The difference signal is processed by a processing 
algorithm, which matricial equation is described in equation 39. The modulating signal is used to 
generate the calculation matrix coefficients after frequency decimation. The initial values is set 
up to avoid overshoot. For small initial values, the overshoot reaches high values, otherwise, for high 
initial values, the setting time becomes very long and the system turn very slow. The ideal values can 
be foun by simulation. The best values of modulation factor (3 m minimizes the scale factor error and 
serves to maintain the modulation range limit into the FSR 3 interval. The output can be sent to an 
adaptive LMS filter algorithm (see reference [10]) to enhance the readout signal or at the digitalized 
input signal, or a predictive Kalman filter (an LMS adaptive filter scheme for the enhancement signal 
configuration is depicted in figure 6). 

4. Simulation Results 

The simulations were made over a FORG wich uses a 10 meter lenght single mode, polarisation 
maintaining, fibre-optic, and 1.446 refraction index and 1.55/xm source lenght. The loop has 10 
centimeter diameter, and the phase modulator is working at the modulation frequency equal 5.186722 
MHz over the modulation factor f3 m ss 1.7278759594 . The sampling frequency is 100 times greather 
than the modulation frequency. These parameters yields a scale factor error is around of 1.18 x 10~ 3 % 
at 20 radians per second is found. All these values are used in this work. The Signal response was 
plotted at 20rad/sec rotation rate within 200 samples (see figure 7), where I CC w(k) counterclockwise 
signal output, I cw (k) clockwise signal output, and signal difference. These signals had been measured 
before pass by decimation frequency process. 

In Figure 8, we can see the signal response at 20rad/sec rotation rate within 200 samples, where 
Iccw{k) counterclockwise signal output, I cw (k) clockwise signal output, and signal difference. These 
signals had been measured after pass by zero order holding and decimation frequency process, and in 
the figure 9, we have the signal response at 20rad/sec rotation rate, the signal output at 200 samples, 
and the signal's FFT and phase response. 

Reverting the rotation sense, note a signal response invertion (at -20rad/sec rotation rate: signal 
output at 200 samples, FFT and phase response). Note the phase inversion in relation of the previous 
figure. At -2rad/sec (depicted in figure 11), we get a signal response (showing signal output, FFT 



and phase response at 200 samples only), we can note that the FFT response varies proportionally 
(compare with the previous figure). It shows that the alghoritm is working properly and the equations 
shown in section 2 may appear valid. 

But we can note the residual ripple due to the higher order expansion error approximation 
too (see figure 12). These riple can be easily minimized by low pass digital filtering. This fluctuation 
introduces a coloured noise, causing a bias drift at a mean value, masking the signal. 



5. Conclusions 

The use of an optical switching and digital frequency decimation can simplify the the optical 
and electronic hardware, reducing the photodiode offset, getting a signal directly proportional to the 
rotation rate with sense information. It can do the digital spectrometry processing method as an 
useful way to improve a low cost resonant gyroscope performance using less fiber. In this technique, 
the FORG performance is dependent of modulation factor (3 m , as well as the decimation factor M, 
that smoot the response curve and reduces the processor's calculations cost. 
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